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, Cayley tree. In the limit of zero field, the system displays glassy behavior below a crossover 

I temperature, Tg, that scales inversely with the logarithm of the number of generations; 
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Abstract 

We present a detailed study of the nearest-neighbor ferromagnetic Ising model on a 



thus Tg is inversely proportional to the logarithm of the logarithm of the number of 
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sites. Non-Gaussian magnetization distributions are observed for T < Tg, reminiscent of 

I that associated with the central spin of the Edwards- Anderson model on the same tree; 

C ' furthermore a dynamical study indicates metastability, long relaxation times and ageing 

. consistent with the development of glassy behavior for a finite but macroscopic number 

SIh i of sites. 
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1 Introduction 



Recursive structures like the Bethe lattice and the Cayley tree provide a pedagogical enviroment for 
the study of physical problems; in this setting they can be treated with a direct analytic approach 
without resorting to approximate methods. The Bethe lattice, an infinite Cayley tree, is a connected 
dendritic structure with constant coordination, z, and no loops, as displayed for 2; = 3 in Figure |l|. 
Strictly speaking it is a pseudo-lattice since it cannot be embedded in any real finite-dimensional 
lattice; indeed it is often regarded as an infinite-dimensional structure since the number of sites 
accessible in N steps from a given site (~ N"^ for a d-dimensional lattice) increases exponentially 
with A^. Thus the Bethe lattice provides a setting where mean-field treatments can become exact. 
This property was first discussed by Domb who showed that the Bethe-Peierls (BP) approximation 
to the nearest-neighbor (nn) ferromagnetic (FM) Ising problem, with H = —JJ2{ij)'^i'^j where 
J > 0, (Tj = ±1 and (ij) indicates a nn sum, is exact on this structure; IH its solution is identical 
to that of the infinite-range FM Ising model. Similarly Thouless, Anderson and Palmer studied 
the infinite-range Sherrington-Kirkpatrick (SK) model of spin glasses on the Bethe lattice using a 
mean- field technique; |^] they were able to recover the key results of Sherrington and Kirkpatrick 
(SK) without using the replica method. More recently there have been several studies of the SK 
model on the Bethe lattice, particularly in finite fields. [0, ^, ^ |ll|, |l2|, In general the 

study of a variety of problems on this recursive structure has helped to develop our understanding of 
diverse physical phenomena including self- avoiding polymers, random resistor networks] 17] and 
percolation, [p^] 

Like the Bethe lattice, a Cayley tree is a connected structure with fixed coordination number 
and no loops; however it has a finite number of generations (cf. Figure [l|) and hence sites that 
are dominated by the boundary. More specifically, the total number of sites in a Cayley tree of n 
generations with coordination z is 



N = l + z + ziz-l) + ...z{z - IT-' = (1) 



and the number of surface atoms is 



so that for large n 



Ns = z{z - ir-' (2) 
^ (z-2) 

aT - (731) ^'^) 

in contrast to the situation in "real" lattices ~ N~d^. Thus the "interior" of a Cayley tree, in the 
limit of a large number of generations, contains an arbitrarily small fraction of its total number of sites, 
and the boundary plays a key role in any problem studied on this graph. In particular, the Bethe- 
Peierls transition for the FM Ising model on a Cayley tree occurs only for its central spin; despite its 
finite moment, the total spontaneous magnetization of all the spins remains zero.[|^, |2^, ^ In a 
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nutshell this occurs because, at zero field and low temperatures, very large domains of flipped spins 
can nucleate from the boundaries; the resulting finite-size glassiness is the subject of this Paper. 

The recursive structure of the Cayley tree permits a detailed analysis of the single-site mag- 
netization distribution as a function of field and generation. In doing so, we find that for fields 
h < hco, where the crossover field /ico decreases exponentially with the number of generations, there 
is a temperature-scale Tg below which well-defined, large domains of flipped spins exist. For T <Tg 
the magnetization distribution becomes non-Gaussian, reminiscent of that associated with the cen- 
tral spin of the magnetized spin glass phase of the itJ model on the Cayley tree. The cross-over 
temperature, Tg, scales inversely with the logarithm of the number of generations of the Cayley tree 
so that the "finite-size" glassiness persists to very large system sizes for a macroscopic number of 
sites. 

We therefore have a short-ranged periodic spin model that has a "glass" cross-over temperature 
that decreases very slowly with increasing system size; more specifically it is inversely proportional 
to the logarithm of the logarithm of the number of sites. We characterize it using a combination 
of analytic and numerical techniques, always retaining open boundary conditions. First we study 
the magnetization for different thermodynamic limits emphasizing the crucial role of the ratio of 
surface/bulk sites as n ^ oo. We recover the Bethe-Peierls result if this ratio goes to zero; otherwise, 
for vanishing applied fields, there is a cross-over to a glassy phase characterized by well-defined 
clusters of flipped spins. We find in that for h < h^o, the single-site magnetization distribution 
becomes non-Gaussian for T < Tg similar to that of spin glass models residing on the same structure; 
however it recovers its Gaussian character with increasing field. For T < Tg the largest barriers 
associated with developing broken bonds in these domains scale with the number of cluster sites; 
we thus refer to this low-temperature state as a finite-size glass. A dynamical study of this system, 
performed numerically, indicates the presence of metastable states and long relaxation times at low 
temperatures. The autocorrelations for T < Tg are determined after a waiting time, and indicate 
ageing effects; the variation of x'l and x's with temperature also agree with the presence of glassiness. 
As expected, the Edwards- Anderson susceptibility of the entire tree has a maximum which develops 
slowly with system-size; no divergence is observed. Finally we find that for this system the density 
of Lee- Yang zeroes in the presence of a complex field is very high in the vicinity of the real axis 
below a particular temperature-scale; we identify it as the crossover temperature Tg. We end with a 
summary of our results and plans for future work. 

2 The Different Thermodynamic Limits on the Cayley Tree 

2.1 Warm-Up: The Bethe-Peierls Transition of the Central Spin 

One way to take the thermodynamic limit on an n-generation Cayley tree is to look solely at the 
properties of its central spin, and then to take the limit n — > -|-oo. As first pointed out in |^, the 
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behavior of the central spin is then characteristic of an infinite-dimensional lattice; more specifically 
it displays a mean-field transition. In other words, the Bethe-Peierls approximation becomes exact 
on a finite Cayley tree if and only if one considers solely the properties of its central spin and ignores 
its surface. For recent results in this field, we refer the reader to reference [jl6| In appendix A we give 
the calculation of the recursion relation for the partition function Zn{P, H, Hn) of an n-generation 
tree with coordination z; here /? is the inverse temperature and H and Hn are the magnetic fields 
acting on the spins of generations to n — 1 and n respectively. The resulting recursion relation for 
is 

(^-1)" 

Z„(/?,/i,/i„) = (4cosh2(/3J)cosh2(/?/i„)) ' Zn-i{P,h,h + T.hn), (4) 

where the transformation of the magnetic field is 

1 cosh P{J + h) 
^•^=2^^" cosh/3(J-/.) - 

We consider the special case of a field e applied at the surface of a tree, and ask whether it is amplified 
in the bulk as determined by the recursion relation (^). This condition defines the bulk critical 
temperature f3c by the expression (z — l)tanh/3cJ =1- If /3 < /3c, the magnetization of the central 
spin is zero; however if /3 > /?c, there is a broken symmetry for the central spin in the thermodynamic 
limit. As discussed in Q, the critical behavior of the central spin in the thermodynamic limit is 
identical to that of an infinite-range ferromagnet where [5 = 1/2 and 5 = 3. 

2.2 Beyond the Bethe Peierls Regime: Different Thermodynamic Limits 

We now wish to look at the transition, not only of the central spin, but of the entire tree. First we 
consider only half- space-trees, that is trees such that the coordination associated with the ancestor is 
z — 1 and not z. In what follows, we denote a half-space-tree of n generations as a n-half-space-tree; 
one is pictured on figure |^. We shall label the generations so that the ancestor is at generation n 
and the leaves are at generation 1. We are interested in the magnetic properties of the spins in 
generations n — m to n in the limit n — > +00 where, of course, m is a function of n; more specifically 
we want to classify the different regimes as a function of m(n) in an external uniform field h. In 
order to obtain the magnetization of the generations n — m to n, we must apply a source magnetic 
field A to these generations; then we differentiate the partition function with respect to A in the limit 
A ^ to obtain 

(M(n,m,/i)) = ^^lnZ(n,m,/i,A = 0). (6) 

Details of the calculation of the partition function are given in appendix B, and we obtain 

,Nv x-^ / -xr, i sinh/?/ij cosh/3/ij d/i;,, , , ^, dhn,. 

(M(„,™,.))= E - ')"- ;j ^ ^ = °) + ta..h;9fe„— (A ^ 0). (7) 



i=n—m 
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The last term is the contribution of the ancestor to the average magnetization; hi is the total field 
at generation i which is the sum of the external, the source and the recursive (of (|5|) fields {h,X and 
^jnd-j respectively. 

We do not treat the iteration of /i*"*^ exactly, but approximate it as described below. From the 
iteration (^) we can deduce the shape of hf^i as a function of hi, which is plotted on figure ^. If 
P < (3c, the slope at the origin is less than unity, whereas it is larger than unity for (5 > f3c- Moreover 
for P > Pc there is one non-trivial fixed point, for which hf^^ = hi = h* which depends only on the 
temperature. This behavior suggests that the iteration can be approximated by linearizing 
in the vicinity of hi = and of hi = h*. More precisely, we deduce from (P) that 

dhfli _ sinh/3Jcosh/3J 
dhn cosh^ /3 J + sinh^ phi' 

and define new variables, r]i and r/2, such that 

-^{hi = 0) = (z-l)tanh/3J = l + r?i (9) 

-^(^h, = h*) = l-r/2, (10) 

where r/i G [0, z — 2] and 772 G [0, 1] and we note that r/i and 7^2 depend only on the temperature. We 
then express the complete recursion by 

h^ = {l+rii)h, (11) 

or 

h^iX\ = mh* + {l-m)hi. (12) 
The first linearization corresponds to hi G [0, /ic] and the second one to hi € [h^ h*], where 

he = -^^h*. (13) 

If A = 0, iteration of the total magnetic field hn leads to different results depending on the relative 
magnitude of i (generation) compared with 



nc{h) 



1 + i 



.In (1 + r/i) \r]i + r]2 h 
where [] denotes the integer part; for i < ndh) 



(14) 



/i, = -((l + r?i)^+i-l) (15) 



and if i > ndh). 

hi = h* + - + -h*--]{i- my-'"'. (16) 

There are therefore three regimes which we can study corresponding to the application of large 
(I), intermediate (II) and vanishing (III) fields on the finite Cayley tree; more specifically they 
correspondand to the conditions 0<nc<n — m(I),n — m + l < < n (II) and Uc > n (III) 



respectively where is defined in equation (14). 
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2.2.1 Regime I: < ric < n — m 

We now determine the average magnetization as given by equation (0) in Regime I. Using the ap- 
proximation discussed in the previous section, we write 



^(A = o) = l(i-(i-7?2r"+"+i) 



If < /ij < /ic or i < Tic, we can approximate 

sinh/3/ij cosh I3hi 



cosh PJ + sinh Phi 



aihi 



where 



and if hi > or n > ric, we use 



ai 



cosh pj 



sinh phi cosh phi 
cosh^ pj + sinh^ 



02 (/ij - h*) + 62 



where 



a2 
62 



cosh^ pJ cosh^ + sinh^ pj sinh^ 
(cosh^ /?J + sinh2 /3/i*)2 
sinh/3/i* cosh/3/i* 



cosh^ pJ + sinh^ 



(17) 

(18) 
(19) 

(20) 

(21) 
(22) 



It is now straightforward to insert these expressions into the equation (|^) for the average magnetiza- 
tion. Since we normalize by the number of sites 



l + (z-l) + ... + (z-l)" 



{z - l)™+i - 1 



z-2 

it is reasonable to neglect the contribution of the ancestor. We obtain 



{M{n,m,h)) = — 



{b2 + a2 — )Nm 



(23) 



(24) 



l-((l-r/2)/(^-l))" 

l-((l-772)/(^-l)) 



-(62 + a2-)(l - m) + a2(V - /ioo)(l - r?2)"-'"-"^ 



a2{hn^-hoo){z-lT(l-mr 



^^l l-((l-7?2)V(z-l))"' 
1-((1-^2)V(^-1)) 



where /too = h* + h/r]2. We can now discuss the different thermodynamic limits. First we note that 
if n — > 00, m — > 00 and n — m — > 00, the thermodynamic limit of the normalized magnetization reads 



lim 



{M{n,m,h)) 



62 + 02- 



h 



z - 1 



772/ z-2^r]2 



(25) 
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However if we take a different thermodynamic limit, with n — m ^ b, where 5 is a constant thickness 
boundary, we obtain 



{Min,m,h)) ( h\ z-l (^_i)(^_2)(l-^2)^-"^ 
77 = 02+02 — — a2\n'oo — /In J 7 ^7- 77 r. ^ 



The absolute value of the corrective term decreases as the thickness of the boundary increases since 
II-T72I < 1. 

2.2.2 Regime II: n — m + 1 < < n 

Again we calculate the average magnetization, (^), this time in the regime n — m + l<nc<n where 
we use the relations ( [l8| ) and (^0|) from the previous section. Furthermore we need to approximate 
the expression dhi/dX for A = 0. n — m < i < ric, we have 

^^'i (A = 0) = - ((1 + r?i)-"+-+^ - 1) , (27) 



and if i > 



where 



We have used 



dX rji 



^(A = 0) = l + a(l-r?2r-, (28) 
uA r/2 

a = — (1 + rj^)n,-n+m+i _]__]_ (29) 

--h* + ^ + {i- mr^^ (kA^) -h*- ^] , (30) 

V(A) = KS^) + - ((1 +r?i)"-"+-+i - 1) . (31) 



with 

V(A) = /i„.(0) + 

Vi 

We insert these expressions into and sum the full geometric series. We note that the thermody- 
namic limit can only be taken with n — m^b, where again 6 is a constant. We find that the dominant 
behavior at small magnetic field h depends on the temperature. Let T' be the temperature such as 
tanh/3' = l/\/z — 1. Then if T' < T < Tc the dominant term in the normalized magnetization is 
linear in h with corrections in /i", where a = In (z — l)/ln(l + r/i) — 1; by contrast if T < T' the 
leading term is of order h". We do not write explicitly the corresponding expressions since they are 
tedious and do not contribute further to the present discussion. 

2.2.3 Regime III: ric > n 

In this regime, we approximate the average magnetization by inserting ([T^ ) and (^7|) into (0) which 
yields 

(M(n,m,/i))= (z-ir-Wh^-^^^^ n^+Vi) (33) 



i=n—m 
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After summing the geometrical series, we obtain 



{M{n,m,h)) 
N{m) 



aih 



ryf 1 - (z - 



(z-1) 



1 



z-2 



(l+r?i)((l+r?i)^ + l) 

2-2-7/1 



1 + 

z - 1 



z-l-(l + r?i)2 ^ 
We identify the cross-over field hco as /in-i 



z-1 



he, that is 



(33) 



(34) 



(??l+??2)((l+m)"-l)' 

which is exponentially small in n. If we take the limit of large n and m but finite n — m, we can write 



{M{n,m, hco)) 
N{m) 



h' 



+ 



aim 



(l + Vi)'^' 



z-2){l + m)-- 



1 (l + r,i)((l+??i)'' + l) 



z-2-rii 



1 



z - 1 



1 



z - 1 - (l + r/i)2 \' \ z-1 

In the limit T ^ 0, (1 + f?i)/(z — 1) = tanh /? — > 1, so that the only non-negligible term is 

(l + r?i)-"(l + r?i)2™ 

— = (1 + fyi) 



(35) 



(z-1)^ 



z-1 



Since 



(36) 



(37) 



at small temperature, we obtain the cross-over temperature scale Tg = 2/lnm. In the next section 
we shall rederive this crossover temperature with another method and give its significance. 

3 Magnetization Distribution and Finite-Size Effects 

In this section we study the distribution of the magnetization which can be determined from exact 
recursion relations. We find that finite-size effects are crucial in this analysis, and we recover the 
crossover temperature Tg discussed in the previous section. We also compute the distribution of the 
magnetization in finite field. 

3.1 Magnetization distribution in zero field 

We begin with the case of zero magnetic field. Given z — 1 n-half-space trees with coordination z, 
it is straighforward to obtain a (n + l)-half space tree with the same coordination. One just has to 
add a common ancestor and to link it ancestor to the z — 1 ancestors of each n-tree (see figure |2|) . 
In order to get a full tree, one has to "glue" z half-space-trees instead of z — 1 at the last step. Let 
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Pn{M) be the conditional probabihty for a n-half-space-tree to have a magnetization M, given the 
spin of the ancestor is a. Of course, 

E^nW = l (38) 

M 

The recursion relation for P^(M) is 

P^{M) = KM-{Mi + ... + M,^i + a)) 

Nh,...,Mz-i 

E I - ^r''~' n ^n-i(^^o n K^m), (39) 

A,-0 \ / j=l i=fc+l 

where x is the probability for breaking one bond: 

-f3J 



e 



(40) 



e/3J + e-/3J- 

The initialization of the recursion is given by Pf(c7') = 5a,a' ■ This recursion can be performed 
numerically, at least for a small number of generations. The result is plotted on figure ^ for z = 3 
and 10 generations. For a finite size tree and at low temperature, the magnetization distribution 
presents a non gaussian structure, reminiscent of the magnetization distribution of the central spin 
in Bethe lattice spin glasses |10| |14]. Notice that the temperature which controls the departure from 
the Gaussian distribution is lower than the bulk transition temperature. 

In order to determine precisely this temperature, we compute the recursion relations for the 
average magnetization. Using the recursion relations (^), we find 

z-l / 

k 



mUi = 1 + E ( \ - ^y~'~' iHN)n +{z-k- 1)(M)+) (41) 



fc=0 

It is clear that, since Pi*"(<t) = Pf (-o-), for all n P+(M) = p-(-M), so that (M)+ + (M)- = 0. 
Putting this equation into (p|) and using the well-known relations for the sum of binomial series, we 
obtain {M)^_^-^ = p{M)^ + 1 and {M)q = 1 where p is defined as p = (z — 1) tanh (3 J. This recursion 
can be easily solved, and one gets 

im ^ z - 2 - 1 

Nn 1 (Z- 1)"+1 - 1' ^ ^ 



where Nn is the number of sites of a n-half-space-tree, given by (|23|). The structure of the distribution 
of magnetization is non-Gaussian provided 2xn <C 1, that is 

^ ^ " h^ " In (lniV„/ln(z- 1))' '•^^^ 

It is clear that the temperature Tg (which will be identified wit the glass temperature) decreases very 
slowly with the system size. For instance, in limit where ~ 6.02.10^^ and z = 3, Tg = J/4.4. We 



8 



conclude from this analysis that finite-size effects persist in the limit of a macroscopic number of 
sites. For Tg to be drastically reduced, one should consider systems of size exp 6.02.10^^ (!). Since 
the appearance of glassiness is a finite-size effect, we stress that Tg is a cross-over temperature scale 
even for macroscopic systems. 

3.2 Structure of the Magnetization Distribution for T < Tg 

We would like to understand qualitatively the structure of the maxima of the magnetization prob- 
ability distribution below Tg] more specifically we want localize the maxima and to calculate their 
weight. To do this, we use the normalized continuous magnetization variable m = M/Nn S [—1, 1] 
and the associated density Pn{m) = NnP^i^)- The recursion relations for pn{m) are derived in 
a straightforward fashion from those for P^{M) that are shown in (|39|). Since this relation is a 
convolution, we write the recursion in terms of the Fourier transform Pn{k) of pnirn): 

p^{k) = / e''''^ pn{m)dm (44) 



We obtain the recursion of the p's: 

~Pn+l{k) = [xpn{-^!^) + {l-x)p{-^)^ (45) 

Po{k) = e^' (46) 

We now study the special case z = 3 where the formulae are simpler. It is easy to check by recursion 
that 

n-l . , s 

Pn{k) = (1 - 2x(2" - + 2x ^ 2" exp i(l - —)k] + 0{x^). (47) 

As we shall see later, this low-temperature expansion is meaningful below Tg even in the presence of 
a finite density of kinks. The expansion ( ^7| ) tells us that pn(m) has peaks for = 1 — 1/2", where 
a G (0, n — 1). Moreover, we find that the weight of the peak a + 1 is twice that of a. Inspection 
of figure ^ indicates that this prediction is correct, at least in the region where the overlap between 
the peaks is small. 

These results can also be interpreted in the following fashion. The expansion ( p7[ ) at order x 
means that the magnetization density is calculated at the order of one kink. It is clear that a single 
kink at generation n — a leads to a magnetization 1 — 1/2" and that the number of choices to put 
a kink at generation n — a — 1 is twice the number of choices to put a kink at generation n — a, 



which is the content of equation (47). What is striking is that below Tg, this one kink picture is 
valid, even though we deal with a finite density of kinks x. This means that Tg is the temperature 
below which the kinks are rarely nested. In order to check this assertion, we determine the condition 
for when the kinks to "induce" well-defined domains of flipped spins; this criterion will determine 



the validity of the expansion (|47|). The overlap between the domains induced by the kinks is small 
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provided xNn{S)n < Nn- In this expression, {S)n is the average size of a domain of flipped spins 
induced by a single kink. The number of descendants of a kink at level n — p for a n-half-space-tree 
is ^p.n = 1 + 2 + ... + 2"-P = 2"-P+^ - 1. The average over p of Sp^n is 

= ^g_',r = - 1 - n - 1. (48) 

The condition for "non-overlap" is xn < 1 which is just T < Tg. We conclude that for T < Tg, the 
number of kinks is small enough for the system to develop well-defined domains of flipped spins. In 
this temperature regime, the excitations of the spin system are kinks which are local in the bond 
variables but are highly non-local in terms of the spins, in terms of bounds, but highly non-local in 
terms of spins. 

3.3 Magnetization Distribution in a Magnetic Field 

The magnetization distribution in a magnetic field can be computed using 



(m) - f49) 



The conditional magnetization distribution (the central spin being parallel to the field) is found 
to converge rapidly to a gaussian as the field increases. Since large domains of flipped spins do not 
survive in a magnetic field, the cross-over field is expected to decrease drastically with the system size. 



Indeed, we have shown previously that the cross-over field hco of equation (34) decreases exponentially 
with the number of generations. The distribution of magnetizations is plotted on figure ||, where the 
central spin is taken to be parallel to the field. 



4 Barrier Structure 



We have shown that the excitations below Tg are broken bonds. In order to characterize the dynamics, 
we calculate the barriers associated with these excitations. The energy barriers at zero temperature 
of a half-space-tree are defined as follows: one starts with a configuration where all the spins are up 
and then one considers single-spin-flip paths from the initial configuration to a final one where all 
the spins are reversed. To each such single-spin-flip path, we associate the maximal energy reached 
during the "passage" from the initial to the final configuration where we take the energy of the 
former to be zero. Then the barrier is defined as the minimum over all the paths of the maximum 
energy of one path. Typically, the Monte Carlo algorithm samples all the paths in an ergodic way; 
by contrast the Swenson algorithm does not generate paths with respect to the single spin flip. In 
Appendix B we give details of the calculation of the barriers associated with the n-half-space-tree. 
Here we calculate the number of states with a given barrier at a given temperature below Tg, 
for a full n-tree. Following our previously established convention, the leaves and the center reside at 
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the first and the n-th generation respectively. We note there exist n* kinks at a given temperature 
where n* = Nx. The barrier for a configuration of n* kinks is assumed to be only a function of 
the generation a of the kink which is the closest to the origin. In order to calculate the number of 
states with an energy barrier E^, we have to enumerate all the configurations with no kink between 
generation a + 1 and n, kinks at generation a and n* — Ua kinks between generations a — 1 and 
1. We call the number of such configurations of kinks g{a). We make the approximation that the 
energy barrier of all these configurations is so that its life-time is, according the the Arrhenius 
law, 

Ta = To exp {xpEi), (50) 
where A is a constant. Such a configuration of kinks is displayed on figure H Clearly, we have 

n* 

^(a) = 9a{na), (51) 

na=l 

with 

We can calculate the sum and obtain 

g{a) = ~ l)'^"""')^ _ l^^2 ((- - 1)""'^- - 1)"-")^ ^ (53) 

and we get the probability P{a) for the system to be in a valley with a barrier Ea. 

-1 r 



iV(l -{z- l)-")\ _ /iV(l -{z- 1)-"+!)' 
Nx I ~ \ Nx 



(54) 



where we have normalized by the total number of accessible states at a given temperature on a 
n— tree. We have also approximated the number of sites as follows: 

N = ^ ~ ^ . (55) 

z-2 z-2 ^ ' 

Assuming both that x <C 1, and that {z — 1)~" ^ 1, and using Stirling's formula, we get 

^(«) - exp [-j^) - exp (- (7^^) • (56) 
If Nz/{z - lY <C 1, we obtain 

Pia) ~ Nx ,^~\ . (57) 

As expected, P{ot) decreases as a function of a, which indicates that the long-lived states are less 
numerous than their short-lived counterparts. 
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5 Glauber Dynamics 

We now study the Glauber dynamics of the spin system on the Cayley tree. 



5.1 The Glauber Matrix 

We begin with a general discussion of Glauber dynamics. Let p{{a},t) be the probability for the 
system to be in a state {a} at time t. Following Glauber |23|, we define Wi{{a}) as the probability 
per unit time that the spin i flips from cij to — Cj, while the others remain fixed. The master equation 
is 



/ N \ N 

[^Wi{{a}) 1 p{{a},t) + ^Wi{{ai,..., -ai, ...,aN})pi{(7i, -ai, ctat}, t). (58) 

\i=l / i=l 



Since we require the Boltzman distribution to be a fixed point, the coefficients Wi{{a}) are of the 
form 

Wi{{a}) = ^{l- a, tanhip J ^ a,)), (59) 

jevii) 

where V{i) is the set of neighbors of the site i. If one denotes by p{t) the 2-^^ vector of the p{{a},t), 
equation (|5^) can be written as 

|p(i) = G.p(i), (60) 

where G is the Glauber matrix. We first show some properties of the matrix G. Since the Boltzman 
distribution is a steady state of the dynamics, its corresponding eigenvalue is zero whatever the 
temperature. Even though it is not symmetric, the matrix G can be diagonalized and its eigenvalues 
are real; we give the proof of this statement here. The Glauber matrix satisfies detailed balance 
which means that Ga^^pjj'^ = G/j^^pa ^ where p^''-' is the Boltzman distribution. As a consequence 

Let us define a matrix M such that 

M^, = {pi^f'G^,{pfy\ (62) 
Then M is symmetric. Let p be a right eigenvector of the Glauber matrix. Then 

Ga/3P/3 = Ap/3 (63) 

is equivalent to 

E^a,{pfy'^\, = X{p^^^y'^"pa, (64) 
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so that I Pa ) Pa is an eigenvector of M. We conclude that G is diagonahzable, and that all of 



its eigenvalues arc real. 

The spectrum in the infinite temperature limit can be understood in the following manner. If we 
define a state vector by 

IV')=E/((W)l^i)®-®l^^)' (65) 

then its dynamics are 

7 AT 1 A'' 

^l^) = -ylV') + 2E<lV'), (66) 

i=l 

SO that the eigenvalues of the Glauber matrix at infinite temperature are of the form 

^ = -y + ^E/^- (67) 

where Hi = ±1. The spectrum in the infinite temperature limit is composed of levels at integer values 
between —N and 0, with a degeneracy given by the binomial coefficients. 

For bipartite lattices, such as the square lattice of the Caley tree, the spectrum of the matrix G 
is symmetric; more specifically if A belongs to the spectrum, then — — A is an eigenvalue too. We 
give the proof of this statement now. Let Xja} be an eigenvector of M, with an eigenvalue A: 

N ^ I 

XX{a} = - ^ -(1 - tanh {pJhi)X{a} + ^ ^ hJjh-^^'^^' ^^^^ 

where hi is defined by 

hi= Yl "j- (69) 

jev{t) 

Let Yjcr} be defined as 

Y{a} = (-1)-Mx{a}, (70) 

where ^{cr} is the number of up spins in the configuration {a}, {a} is deduced from {a} by flipping 
the spins of one of the two sublattices. Then, 

(MY)M = -g-(l-a,tanh(/3g/.,))(-irWx{a} + E ^ 2cosh(^J^,) -a,, 

^1 _ 1 



= (-l)^M 
= -iN + X)i-iy^''^X{a} = -iN + X)Y{a}. 



g - (1 + a, tanh im^)) X{a} - g ^^^^^^X{ai, -a„ a^} 



Given an eigenvector X for the eigenvalue A, we have constructedd an eigenvector Y for the eigenvalue 
-AT- A. 
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The difference between (^0|) and the Schrodinger equation is that quantum mechanics preserves 
the scalar product which results in Hermitian Hamiltonians. Furthermore, the physical states "reside" 
in a Hilbert space, and each state of this Hilbert space is physical. In the case of the Glauber matrix 
there is no such vectorial space; more specifically, the sum of two probability distributions is not a 
probability distribution. However some quantities are conserved by the dynamics. It is easy to show 
that the eigenvectors of G for the non-zero eigenvalues have the property that 

E P{^} = 0- (73) 

This is a simple consequence of the fact that the Glauber matrix preserves the quantity 

Ep{^}- (74) 

{^} 

5.2 Spectrum of the Glauber Matrix 

We now discuss the spectrum of the Glauber dynamics as a function of temperature. We have 
performed a series of numerical diagonalizations on small clusters, and then speculate on the behavior 
for larger systems specifically with respect to the statistics of their eigenvalues. Since this aspect of 
the Glauber dynamics has not studied in much detail before, we examine the Ising model on one- 
and two- dimensional lattices with this approach before ending with the Cayley tree; in particular, 
we examine their relaxation spectra as a function of temperature and attempt to extrapolate their 
behavior in the thermodynamic limit. 

We begin with the Ising chain, where we have determined the spectrum as a function of tempera- 
ture for a six-site system with open boundary conditions; the result is displayed in figure 0. We also 
treated the case of a 10 sites ring. Of course, we did not follow the 1024 levels as a function of tem- 
perature, but computed the eigenvalues for a given temperature. The spectrum is plotted on figure 
^. From these results, we can speculate on its behavior for large systems. The spectrum is completly 
integrable in the sense that no level repulsion occurs in the evolution of the eigenvalues as a function 
of temperature. It is clear from figure ^ that the evolution of 1/rj is monotonous as a function of 
f3 and that no repulsion is observed. However, the distribution of eigenvalues is non uniform, since 
a curvature can be seen on figure ^. It is not surprising that the spectrum reveals an underlying 
integrable dynamics since the dynamics was indeed integrated by Glauber in [^]. We also note the 
existence of a gap at non-zero temperature, which vanishes in the limit of zero temperature; here 
the eigenvalue is doubly degenerate, corresponding to the existence of a symmetry-breaking phase at 
T = 0. 

We have also studied small clusters (of size 3x3) of the two-dimensional Ising model; the resulting 
spectrum in the low-temperature phase is plotted on figure ^ The obvious difference with the results 
for the one-dimensional chain is the presence of clusters of relaxation times in the low-temperature 
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spectrum. The evolution of the eigenvalues as a function of the inverse temperature is plotted for a 
small 3x2 cluster on figure Level repulsion is visible in the vicinity of the transition temperature. 
For an infinite size system, one expects the presence of a gap between zero and the first negative 
eigenvalue in the spectrum for T > T^, it should vanish at the transition temperature, leading to 
a doubly-degenerate zero eigenvalue in the interval [0, Tc] which corresponds to the existence of a 
symmetry-breaking stable phase. Thinking of the spectrum as a function of (3, one expects the 
spectrum to display chaotic behavior in the vicinity of the transition temperature, leading to a 
rearrangement of the inverse relaxation times at low temperature and the presence of clusters of 
relaxation times in the zero-temperature limit. However, it is not clear whether the level-spacing 
statistics should be Gaussian Orthogonal Ensemble (G.O.E.) in character in an appropriate sector; 
one might be left with a mixture of G.O.E. and Poisson statistics. A better understanding of these 
level-spacing statistics would involve the search of the symetry sectors using group theory which will 
be pursued in a future project. 

In our study of the Cayley tree, we determined the spectrum of eigenvalues as a function of 
temperature for z = 3 nad n = 1, which corresponds to four sites. The resulting "spaghetti" are 
plotted in figure Though the system size is very small, one can still see repulsion of the eigenvalues 
in the vicinity of the bulk transition temperature. We determined the spectrum of the Glauber matrix 



at low temperatures for n = 2, resulting in 1024 eigenmodes; the eigenvalues are plotted on figure 12 
The spectrum displays clusters but, unlike the case of the two-dimensional Ising model, these clusters 
are equidistant and occur at integer values. Furthermore the A = cluster contains 8 eigenvalues 
whereas it contained just two for the 2D Ising model. For the case of = 10 and /? = 2, shown in 
figure |l^, we find that n* ~ 0.18 < 1. In the previous section, we conclude from equation(^) that 
P{a) is independent of the temperature; this result is confirmed by the spectrum of the Glauber 
matrix since the A = and the A = — 1 clusters are separated by a gap. 

5.3 Determination of the Glauber Dynamics 

We now diagonalize the Glauber matrix in order to determine the dynamics of the model in an 
explicit fashion. Let Uq be the eigenvectors of G: G.Uq = AqUq , and let P be the passage matrix 
from the natural basis of pure states to the basis Uq,: 

Ua=^'Paf3ep. (75) 

We look for the temporal evolution of the states e^. At the initial time (t=0) 

Pa(0) = e„ = ^P;iu^. (76) 

and at t > 0, the state is a mixture of pure states and is given by 

p„(t) = ^p-ie^^*^P^,e,. (77) 

13 7 
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We can easily compute the magnetization of Pa{t)- We applied this procedure to the case of a tree 



with one generation. The evolution of the magnetization of the 16 pure states is plotted on figure |13 
at low temperature. Of course, on very long time scales, the magnetization of all the states relaxes to 
zero, due to the fact that the eigenvalues associated to the symmetry breaking state are not strictly 
zero. The evolution of the pure states indicate the existence of metastable states, which are the 
precursors of the metastable states present for larger values of the number of generations. 

5.4 The Factorization Approximation 

We recover the bulk critical temperature by considering the asymptotics of a simplified Glauber 
dynamics. Glauber has shown |^] that one can replace the 2^ variables of the linear dynamics by 
a hierarchy of N, nonlinear, coupled equations for the correlations functions. This procedure is very 
similar to the transformation of the Liouville equation into the B.B.G.K.Y. hierarchy in the kinetic 
theory. The first equation of the hierarchy is 

^(?,(t) = -g,(t) + (tanh(/3J Yl ^j))- (78) 

j£V{i) 

In this expression, q = (a) and V{i) is the set of neighbors of site i. In the case of a one dimensional 
chain, one can use the fact that 

(tanh/?J((Ti_i + (Tj+i)) = ^tanh(/3J)(gi_i + Qi+i), (79) 

and one gets a closed equation for the 1-point correlation functions. It is also clear that the whole 
hierarchy decouples, and that one can use this decoupling to integrate the dynamics. In the case of a 
z = 3 tree, one has to take into account the fact that the sites inside the tree have three neighbors;by 
contrast the leaves have one neighbor. For this coordination, 

tanh/3J((Ti + a2 + 0-3) = a{ai + a2 + o^) + 7(0-1 + (T2 + 0-3)^, (80) 

where the coefficients a and 7 are determined by 

a = ^(27tanh/3J-tanh3/3J) (81) 

J = ^(tanh3/3J-3tanh/3J). (82) 

We can thus obtain the first equation of the hierarchy in the case of the z = 3 tree. For the sites 
with three neighbours, 

jQi = -q, + {a + 77) E Qj + 67( J] ^1) (83) 

jev{i) jev{i) 

For the leaves of the tree 

-j^Qi = -Qi + Qj tanh /? J, (84) 
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where j is the neighbour of i. The factorization approximation consists in decouphng the third order 
correlations into 

( n = n ^j- (85) 

j&V{i) j&Vii) 

This approximation leads to the bulk behavior in the high temperature phase and in the vicinity of 
the transition. We start from a configuration of spins such as qi{0) = qj{0) if the sites i and j belong 
to the same generation. Then, for t > 0, qi{t) = qj{t) if we work with the factorized dynamics. The 
factorized dynamics depends only on n variables, one per generation, and is 

^ = -qn + 3(a + 77)g„_i + G-fql_i 

^ = -9i + (a + 77)(2gi-i +gi+i) + 67gf_igi+i (86) 

^ = -gi + g2tanh/3J, 
at 

where i runs from 2 to n — 1. The equilibrium properties are calculated from the asymptotic values of 
the dynamics. We note that the factorized dynamics possesses a non-trivial fixed point for a finite size 
system, whereas the complete dynamics possesses only the Boltzman distribution as a fixed point. 
The asymptotics of the factorized dynamics is found to reproduce quite well the bulk properties of 



the tree. On figure 14, g„(-|-oo) is plotted as a function of the inverse temperature. This curve is in 
agreement with the fact that the central spin exhibits a mean field like transition at f3c — 0.54. We 
also plotted the asymptotic magnetization of the m generations which are the closest to the central 
spin. It is clear that the predictions of the factorized dynamics are qualitatively wrong as soon as one 
goes out from the center. For instance, it is clear that the entire tree does not develop a transition 
at P = (3c- We note that, in equations (|86D, the transition temperature is determined entirely by the 
stability of the zero fixed point of the linear problem. Below Tc, the largest eigenvalue of the linear 
problem is positive, and negative above T^. As in the Ginzburg-Landau theory, the nonlinear terms 
are responsible for the maximum bound on the dynamical variables in the low- temperature phase. 

6 Monte Carlo Dynamics 
6.1 Relaxation of a Single kKink 

We begin by considering the case of a single kink and look for the relaxation of this excitation at low 
temperatures. Let's call K the set of descendants of the kink. At time t = 0, the configuration of 
kinks is such as ai = —1 is i ^ K and = 1 \ii ^ K . We follow the magnetization of the spins of K 



as a function of time, for various size of K. The result is plotted on figure |15|. The rapid relaxation 
at small times is attributed to the fact that the initial state is not thermalized; the thermalization 
occurs at small time scales compared to the collective processes of crossing the barrier. We define 
the typical relaxation time r as the time-scale associated with the vanishing of the magnetization. 
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In Figure 16 we have plotted the logarithm of this relaxation time, In r, as a function of the number 



of generations of K {uk)- The points are approximately aligned, indicating that the Arrhenius law 

r ~ To exp ( /?A^--^n ) (87) 



2 

is well-satisfied. The case A = 1 corresponds to the relaxation of a n-half-space tree where the 
ancestor is free. However, in our case, the ancestor is not free since the domain K is connected to the 
remaining up spins. The barrier height is essentially the same as in the case A = 1, but the number 
of paths to reverse the magnetization is changed. 

6.2 The A.C. Susceptibility 

Experimentally, the a.c. nonlinear susceptibility has been a very useful probe of glassiness,p^] and 
we can apply the same techniques numerically to study the tree problem that we are discussing here. 
In parallel with what is done experimentally, a small a.c. external field H[t) = hsiuLot is applied to 
the spin system, and we measure the magnetization response, expanded into its Fourier components 

M{uj, t) = J2^'k sin kut + 6'^ cos kujt, (88) 

fc>0 

which has only odd harmonics. In practice, we let the system relax for four periods and we take 
measurements only during the fifth one. The in-phase susceptibilities are related to the in-phase 



Fourier coefficients as |24| 



O'l = x[h + lx'sh'+'-x',h' + f^xW + - (89) 
1 91 

e's = jx'sh' + ^x'nh' + '^^X'rh' + ... (90) 

O'n = ^X,h' + ^xW + - (91) 

O7 = ^xW + -- (92) 

Since all the nonlinear susceptibilities are divergent at the critical temperature of a spin glass, we 
must include all the measurable higher order harmonics in order to properly include their influence 



on the lower ones. [ 24 1 



The susceptibilities Xi ^-nd x'3 plotted as a function of temperature in Figures 17 and 18 
respectively for a number of frequencies. Both Xi a-nd x'3 display maxima, where that of the latter 
is more pronounced; in a "real" spin glass in the thermodynamic limit one would expect Xi to have 
a maximum, and X3 to diverge. We note that x'3 is negative near its maximum, which agrees with 
experiments on spin glasses [^] . The frequency dependence of the position of the maximum in x'l and 
x'3 is also consistent with experiment, ||2^ since the temperature for which x'l and x'3 maximum 
increases with frequency. Thus this study strongly suggests that the Ising model on a finite Cayley 
tree exhibits glassy-like behavior at low temperatures; we shall be more specific when we analyze the 
Edwards- Anderson susceptibility of this model. 
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6.3 Autocorrelation Functions and Ageing 

A good test used to determine the presence of glassiness in a given model involves the computation 



of the autocorrelation functions |25] to determine whether they exhibit ageing. The autocorrelation 
functions are 

1 ^ 

C{t,ty,) = — ^((Ji(t + t^)c7i(t^)) - (cTi(t + t^))(cJi(t^)), (93) 
i=l 

where the sample is rapidly quenched below Tg from a disordered high temperature state. The Monte 
Carlo dynamics runs from times i = to t = the autocorrelations are measured at t = tw where 
an average is taken over the initial configurations. 

The autocorrelations in the high-temperature phase are plotted on figure They decrease 
rapidly with the time t and are independent of the waiting time, as expected; this is simply a check 
of our code. Figure ^ represents the autocorrelation functions below Tg. The aucorrrelations increase 
with increasing t^, indicating the presence of ageing. Such behavior has been observed in a wide 
class of glassy models, including the fully frustated hyper cubic model [p^j . 

7 The Edwards- Anderson Order Parameter and Susceptibility 

The recursive structure of the Cayley tree permits us to compute the Edwards-Anderson order 
parameter and susceptibility at all temperatures. This calculation has already been done for the it J 



model with uncorrelated boundary conditions in [10] |14]. These authors found an Almeida-Thouless 



line for their spin glass model, and we wish to use the same technique to analyse the model discussed 
here. In order to compute the Edwards-Anderson order parameter and susceptibility, we consider 
one tree plus one replica, with an intereplica coupling R. If {a} and {a'} are the spin configurations 
of the two replicas, the Hamiltonian reads 

N N 
(Lj) 1=1 i=l 

As per usual with free boundary trees, we have to distinguish between the properties of the central 
spin and the whole tree; in what follows we discuss both cases below. 

7.1 The whole tree 

The partition function can be calculated by the construction of figure even in the presence of the 
intereplica coupling. We take derivatives of the partition function with respect to the intereplica 
coupling, R, resulting in the Edwards- Anderson order parameter 

1 ^ Id 
,,, = -i^^..,.)^-—,nZ(R=0) (95) 
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and the Edwards- Anderson susceptibility: 

1 / ^ ^ \ ^ 8 



(n) 

We now let Z^J be the conditional partition function of a n-half space tree with respect to the 
ancestor's spins a and a'. It is clear that 

^(n+l) ^ e/3Rg2/3//|^g2/3J^W^^W +g-2/3J^(n)^^-l (97) 

^ e/3V'3^(e-2/3^zS + zi"2 + zi"| + e2/^^zi"2)'"\ (100) 

The initial conditions of the recursion are Z^_^_ = e^^e^^^ , — — and Z^_ = 

g/3_Rg-2/3_ff _ For all n, we can show that = zi"^ and we then define zj"^ = Zq"'* = Z^_^l_ = 

Z^^ and Z"]^ = Using this notation, the recursion relations become 

^ g-/3R |^^{n) ^ 3 cosh /? JZj") + Zil^ (102) 

^(-+1) = g/3(i?-2//) (e-2/3J^(n) ^ 24") + e^/^^zi"!)) . (103) 
which can be used to obtain qea and xea in a straight-forward fashion . qea and xsA are plotted 



in figures 21 and 22 respectively for n = 10 and n = 80. In both cases, the curves corresponding to 
these two system sizes are similar, even though the number of sites is small for n = 10 (2047 sites) 
and macroscopic for n = 80 (4.02 moles of sites). We believe that this behavior is related to the 
very slow variation of the glass temperature Tg with the system size. In figure ^ we see clearly that 
Xe.a. does not diverge in the thermodynamic limit but rather exhibits a maximum, characteristic of 
a finite size effect, even in the macroscopic regime. 

7.2 The central spin 



Following the authors of |1C] [14|, we note that 

Zin) ^ zf") + 2Z(") + (104) 
= ^(Z;")-24")+zi")). (105) 

The Edwards- Anderson order parameter of the central spin is qe.a. = Q^"'\ and the Edwards- 
Anderson susceptibility is 
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where the tilde symbol denotes quantities with respect to the central spin. The "central spin" 



Edwards- Anderson susceptibility xe.a. is plotted in figure ^ for various system sizes. Though we 
have seen that finite-size effects are negligeable for the whole tree, they become crucial for the central 
spin. More specifically, the maximum values of the Edwards Anderson susceptibility of the central 
spin does not increase (see figure ^4|) as a function of coordination z. As the system size becomes 
macroscopic, the glass temperature Tg depends more on the coordination than on the number of 



generation; in particular, it increases with increasing coordination consistent with equation (|43|). 
7.3 Central spin in a magnetic field 

In the Sherrington-Kirkpatrick model, glassy behavior persists even in the presence of a magnetic 
field; the Edwards-Anderson susceptibility diverges with an exponent 7 = 1 through the Almeida- 
Thouless line. We have studied the behavior of the Ising n.n. ferromagnet (fm) on a Cayley tree 
ferromagnetic model with free boundary conditions with an applied magnetic field. Specifically we 
studied the variation of the Edwards- Anderson susceptibility in a field as a function of system size. 
As displayed in figure |2^, the maximum of the Edwards- Anderson susceptibility as a function of 



temperature depends on the size. In figure [S^, we showthe locus of the maxima of xe.a. for different 
magnetic fields as a function of temperature. The Edwards- Anderson susceptibility decreases strongly 
with increasing magnetic field, indicating the vanishing of glassy behavior consistent with expectation. 

8 Lee and Yang zeros 

The partition function of the Ising fm on the Cayley tree in the absence of a magnetic field is equal 
to that of the linear chain, so that the Lee and Yang zeros in the plane of exp2/3 have the identical 
structure in both models. Clearly this will not be the case for finite applied field. We can introduce 
a complex magnetic field and determine the zeros of the partition function in the complex plane 
of exp2/3/i. Since we do not impose a uniform magnetic field, but rather one only on the leaves 
of the tree, we can compute the zeros using the recursion relations (|5|), which involve only second 
order equations. Since the magnetic field is non-uniform, the underlying assumptions of the circle 
theorem [^] are not fulfilled. The zeros are not exactly on the unit circle but areb in its vicinity. 
More precisely, numerical work indicates that the zeros approach the unit circle with increasing 



generation. The zeros are plotted on figure ^ in the cases j3 < (3c, P = f3c and (3 > f3c- In figure ^ 
we exhibit the ratio of the density of zeros in the vicinity of exp 2f3h = 1 over the averaged density on 
the circle. The ratio is zero for temperatures larger than the bulk transition temperature, increases 
to a small positive value below the bulk transition temperature, and eventually becomes larger than 
unity below a certain temperature. In the case of a Weiss model, the density of zeros g(9) is given 



by ^ 



9(9) = ^{1 -2r cos (107) 
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where 

r = \J A + r'^ — Ar cos (/9exp ^—2-^(1 — r cos Lp)^ (108) 

?c / rsmp \ 

= — 2— r sin G9 + (/9 + arctan . (109) 

T \2 — r sm cp ) 

The ratio of the density of zeros at the point exp 2/?/i = 1 over the average density of zeros can be 
determined numerically from the sysetem of equations above, and is plotted on the inset of figure 
(^8|), and is always less than unity. By contrast on the tree, the density of zeros in the vicinity 
of the real axis is thus abnormally high below a certain temperature, which we identify with the 
temperature scale Tg. 



9 Discussion 



In summary, we have studied the static and the dynamical properties of the nearest-neighbor Ising 
model on a Cayley tree. At zero field, we find that this system displays glassy behavior below a 
size-dependent temperature that scales inversely with the logarithm of the number of generations; 
thus its glassy behavior persists for a finite but macroscopic number of sites. Because the ratio of 
the number of surface to bulk sites, and the strength of the external field, /i, play a key role 
in the physical behavior of the resident spin system, the different thermodynamic limits associated 
with the values of ^ and h are characterized; the cross-over temperature, T^, is associated with 
fixed ^ in the limit of vanishing applied field. Physically well-defined large domains of flipped 
spins develop at Tg\ at this temperature the probability of nested spin clusters is small. The largest 
energy barriers associated with overturning these domains is determined to scale logarithmically with 
the number of sites at zero temperature, a result that should be valid at finite, low temperatures 
if overlap between spin clusters does not occur. A dynamical study indicates the appearance of 
metastable states and long relaxation times at low temperatures. The autocorrelations are computed 
after a waiting time using Monte Carlo dynamics; they exhibit ageing for T < Tg. The temperature 
variations of the coefficients of Xi ^^id Xs are also determined and they agree with the existence 
of finite-size glassiness. Finally the Edwards- Anderson susceptibility of the entire tree displays a 
maximum (but no divergence) that evolves slowly with increasing system size; that of the central 
spin has much more marked size-dependence. 

We have thus performed a detailed characterization of the low-temperature phase of a short- 
range periodic spin model resident on a Cayley tree. In this particular case, we have found that 
it displays finite-size glassy behavior that remains for a macroscopic number of sites; perhaps it is 
best to characterize this low-temperature phase as a very viscous spin liquid. We note that neither 
intrinsic disorder nor frustration exist due to the initial Ising Hamiltonian; the possibility of many 
low-temperature "cluster" states separated by very high energy barriers is a direct consequence of 
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the unusual topology of the Cayley tree. In many ways we hope that this is a warm-up exercise 
towards the study of spin models on more complicated non-Euclidean lattices, e.g. on a constant 
triangulation associated with a surface of negative curvature, where the intrinsic geometry of the 
host may lead to the possibility of glassiness in the absence of both disorder and frustration. 

10 Appendix A 

We propose another derivation of the magnetization distribution for the tree. Let Z{P, h) the partition 
function of the spin system in the presence of an external field. Then, 

^^^^ = Wj^) ^ " ^ ""^^ (^^"^ '''''' ^ ^ ^ ■ 

Using the Fourier representation of the delta function 



fllll 



we obtain 

where we have used the analytic continuation of the partition function for complex magnetic fields. 
This method is usefull provided one knows how to calculte the partition function, which is certainly 
feasible on a tree. We proceed by decimation, starting from the border of the n-half-space-tree. For 
future purpose, we note /i, the partition function of a n-half-space-tree with a magnetic 

field h acting on the spins of the generations to n — 1, and /i„ on the spins of the generation n. 
Then, 

Z„(/?, h, K) = (4(cosh2 (/5J) + sinh^ {I3K)) ' Z„-i(/3, h,h + T.K), (113) 
where the transformation of the magnetic field reads 

T./. = ^ln^°^^^;-^ + ^;. (114) 
2/3 cosh/?(J-/i) ^ ' 

The last term of the recursion corresponds to the partition function of the ancestor, which is simply 

Zi(/3,/i) = 2cosh(/3r"./i). (115) 

It is straighforward to compute the partition function using these relations and to perform the Fourier 
transform (|112|) in order to obtain the probability distribution of the magnetization. 
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11 Appendix B 



In this appendix we give the value of the energy barrier E{Tn^z) foi" a half-space-tree Tn^z with n 
generations and a coordination number of z for all sites except the root (coordination z — 1) and 
the leaves (coordination 1). We also give the energy barrier E(T*^) for a complete tree T* ^ with 
n generations and a coordination number of z for all sites except the leaves. The derivation of the 



formula is due to A. Sebo and M. Preissmann, and is published in extenso in |29]. Note that the 
same problem arises in the VLSI circuit conception! Generically the problem of finding the lowest 
energy barrier is NP- complete, but the sub-problem of finding the lowest energy barrier of a tree is 
polynomial, and an explicit algorithm is given hereafter. The value of the energy barriers are given 
by: 

c{Tn,z) = r^^^^^l +1 {n,z>3) (116) 

, , , Jn - l)(z - 2), 
c{Tlz) = \- '-^ ^1 + +1 {n,z>2,) 

In the above formula \x\ denotes the lowest of the integers greater than x. The demonstration of 
these formula is constructive. Firstly a lower bound for E is given. Then an algorithm is described 
which produces a labeling of the sites. Flipping the spins in the order of this labeling gives an energy 
barrier exactly equal to the lower bound. The algorithm is recursive. It tries to produce an optimal 
labeling of the sites where the root is labeled before the configuration of highest energy is reached. 
We call strong labeling such a labeling. This extra constraint is useful when one applies z times the 
algorithm on a T^^z to compute E{TnJ^i^z)^ or when one applies the algorithm to Tn^z and to Tn~i^z 
to compute E{T* A strong labeling does not exist when z and n are both odd as shown in [29|. 



Let us now consider on the case of Tn z- Formula 116 means that 



• when z is even the increment in energy when one goes from Tn^z to Tn+i^z is constant and equal 
tof-1 

• when z is odd the increment in energy when one goes from Tn^z to T^+i^^ is alternatively — 1 
and ^ 

We give now the algorithm in the case of even z. The case of odd z is a slightly more complicated, 
but in the same spirit. Consider Tn^z as being made of z — 1 copies of Tn-i^z all of them connected 
to the site 0. Each spin is identified by two numbers k,i with < k < z — 1 and < i < Nn-i 
{Nn is the number of sites of Tn.z)- Let us note vr a strong labeling of Tn-i^z, and uq the root of the 

2 

1. (0,7r(0)), (0,7r(l)),..., (0,7r(Af„_i)) 

2. ... 



(^^2^ — l)th copy of Tn-i^z- The following labeling is a strong labeling of Tn^z- 
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3. (l^-l,vr(iV„)), (£^-l,vr(7V„-l)),..., - 1, vr(no)) 

4. 

5. (^-l,vr(no-l)), - 1, 7r(no - 2)), . . . , (^-1,^0)) 

6. (^,vr(0)), (^,^(1)),..., (^,^(Ar„_,)) 

7. ... 

8. (z-2,7r(0)), (z-2,7r(l)),---, (z - 2, 7r(7V„_i)) 

Note that in step ^, and |^ any admissible permutation can be used instead of the strong labeHng 
vr. It is shown in [^] that the above labeUng is indeed an optimal labeling and it can be used to 
implement a recursive algorithm to find a path between the two ferromagnetic states. 
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Figure captions 



Figure 1: 

(a) A finite Cayley tree and (b) a section of a Bethe lattice with coordination z = 3. 



Figure 2: 

A recursive construction of half-space trees. 



Figure 3: 

/ii4.i(/ii) where hi is the total magnetic field at generation i, and is the iterated field of equation 
(^). The three curves correspond to (a) (3 < (3c, (h) P = (3c and (c) (3 > (3c- 



Figure 4: 

The density probability of the magnetization on a half-space tree where n = 10,2; = 3, the ancestral 
spin is fixed, and the inverse temperature \s (3 = 2>. 



Figure 5: 

The conditional magnetization distribution in the presence of a magnetic field, the field being parallel 
to the central spin. The magnetization distribution is plotted for a z = 3 tree, n = 10 generations, and 
for magnetic fields H = 0,0.001,0.002. The magnetization distribution evolves towards a gaussian 
shape as the magnetic field increases. 



Figure 6: 

A typical configuration of the tree with n = 8 and (3 = 1.5. The vertices with no dots represent up 
spins and the dots represent fiipped spins. Each kink gives rise to a well-defined domain of fiipped 
spins. 



Figure 7: 

Eigenvalues of the Glauber matrix as a function of the inverse temperature for an open six-site Ising 
chain (64 states). The evolution of the energy levels as a function of (3 is monotonous, and the 
eigenvalues are free to cross each other. 



Figure 8: 

Spectrum of the Glauber matrix at a given temperature for a ring of ten sites. The inverse temper- 
ature is /? = 3. No gap is present. 
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Figure 9: 

Spectrum of the Glauber matrix at a given temperature for a 3x3 Ising model with open boundary 
conditions. The inverse temperature is /3 = 3. Some degeneracies appear in the spectrum, and the 
spectrum is symmetric with respect to the A = —N/2 line. 

Figure 10: 

Eigenvalues of the Glauber matrix as a function of the inverse temperature for a 3x2 Ising cluster. 
The number of sites is 6, leading to 64 states. The evolution of the eigenvalues as a function of (3 is 
non monotonous, and avoided crossings are visible. 

Figure 11: 

Eigenvalues of the Glauber matrix as a function of the inverse temperature for a n = 1, 2; = 3 Cay ley 
tree. The number of sites is four, leading to 16 states. The evolution of the energy levels as a function 
of (3 is not monotonous for certain levels, indicating level repulsion. 

Figure 12: 

Spectrum of the Glauber matrix at a given temperature for 2; = 3 Cayley tree with 2 generations. 
The inverse temperature is /3 = 2. The spectrum is symmetric with respect to the A = —N/2 line. 
The clusters collapse at integer values. 

Figure 13: 

Evolution of the magnetization of the 16 pure states of the tree with one generation. The inverse 
temperature is /? = 3. Metastable states appear to be present even for such a small size. The insert 
represents the same curve at short time scales, which shows the transient regime from the natural 
basis of pure states to the metastable states. 

Figure 14: 

Expectation value of the central spin in the simplified dynamics, as a function of the inverse temper- 
ature. The transition temperature is in agreement with the value f3c — 0.54. The shape of the curve 
near the transition is in agreement with the existence of a mean field like transition for the central 
spin. We also plotted the expectation value of the magnetization of the m closest to the origin slices 
of spins, in the factorisation approximation, for m = 3, 4, 5, 6, 7, 8, 9, 10. 

Figure 15: 

Magnetization relaxation of a single kink. Initially, all the spins of a z = 3, n = 10 tree are up, and 
one creates a kink. We call no the number of generations involved in the kink. If no = 1, the kink 
has only one spin, if no = 2, the kink has 7 spins, etc. We follow the magnetization of this domain 
as a function of time. The unit time is one Monte Carlo Step (M.C.S.). One M.C.S. corresponds 
to repeating times the process which consists in choosing one spin at random among the sites, 
and changing or not its direction, according to the Boltzman distribution. The curves are averaged 
over 50 different Monte Carlo runs of the dynamics. 
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Figure 16: 

Logarithm of the relaxation time as a function of the number of generations in the kink. The relax- 
ation time is defined from the cancelation of magnetization on figure [T5|. The points are approximately 
aligned, which is in agreement with the Arrhenius law of equation (p7\j- 



Figure 17: 

x'l susceptibility on the Cayley tree. The tree has 10 generations. The amplitude of the magnetic 
field is 0.1. The curves are averaged over 50 initial configurations of spins, generated at equilibrium. 
One has to find a compromise between the amplitude of the magnetic field and the number of 
configurations to be averaged over, to have a good signal/noise ratio. The curves correspond to a 
period of the magnetic field equal to 500 M.C.S., 1000 M.C.S. and 1500 M.C.S.. 



Figure 18: 

Xs susceptibility on the Cayley tree. The tree has 10 generations. The amplitude of the magnetic 
field is 0.1. The curves are averaged over 50 initial configurations of spins, generated at equilibrium. 
The curves correspond to a period of the magnetic field equal to 500 M.C.S. , 1000 M.C.S. and 1500 
M.C.S. . 
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Figure 19: 

Spin autocorrelation functions of the z = 3 Cayley tree above Tg. The unit time is one M.C.S., the 
inverse temperature is /? = 0.5, and the averages are taken over 100 random initial configurations. 
The autocorrelation functions decrease rapidly with t and are independant on tw 



Figure 20: 

Spin autocorrelation functions of the z = 3 Cayley tree below Tg. The unit time is one M.C.S., the 
inverse temperature is /? = 2., and the averages are taken over 100 random initial configurations. 
The autocorrelation functions depend on the waiting time, and increase with the waiting time, which 
is the signature of glassiness. 



Figure 21: 

Edwards- Anderson order parameter for the whole spin system. The coordination is z = 3. The 
Edwards- Anderson order parameter is plotted as a function of temperature for n = 10 (2047 spins) 
and n = 80 (4.02 moles of spins). The two curves nearly coincide. 



Figure 22: 

Edwards- Anderson susceptibility for the whole spin system. The coordination is z = 3. The Edwards- 
Anderson susceptibility is plotted as a function of temperature for n = 10 (2047 spins) and n = 80 
(4.02 moles of spins). The two curves nearly coincide. 



Figure 23: 

Edwards-Anderson susceptibility for the central spin as a function of temperature for n = 
5, 10, 15, 20, 25, 30, 35, 40 generations. The coordination is z = 3. 



Figure 24: 

Edwards- Anderson susceptibility for the central spin as a function of temperature for z = 3,4, 5. The 
number of generations is n = 20. 



Figure 25: 

Edwards- Anderson susceptibility of the central spin in a magnetic field as a function of temperature, 
for different sizes. The number of generations is n=2, 3, 4, 5, 6, 7, 8, 9, 10, 30. The magnetic field 
is H = 0.05 
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Figure 26: 

Maxima of xea as a function of temperature for different values of the magnetic field. The system 
size is n = 10 generations. The magnetic fields are H = 0.025, H = 0.05, H = 0.1, H = 0.15, H = 0.2. 



Figure 27: 

Zeros of Lee and Yang at different temperatures. The radial coordinate is rescaled in order to allow 
the superposition of different maps of zeros. The external map represents the zeros for 12 generations, 
z = 3 and (3 = 0.3. In this case, the zeros are localized on well defined areas of the complex plane, 
far from the real axis. The intermediate map represents the case 13 = jSc — 0.54. The zeros begin to 
fill the circle, with areas of zero density, espacially at the intersection of the real axis with the unit 
circle. The inner map represent the zeros at low temperature, (3 = 2. In this case, the circle is filled 
with zeros. An analysis of the density of zeros reveals an annormally high density in the vicinity of 
the point h = 0. 



Figure 28: 

Ratio of the density of zeros in the vicinity of the point exp 2/3/i = over the average density of 
zeros in the case of the z = 3 tree with 17 generations. The inset represents the same quantity for 
a Weiss model with = 1. In the case of the tree, the ration is greater than one below a certain 
temperature, whereas it is always less than unity for the Weiss model. 
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